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ABSTRACT 


A unidirectional fiberglass epoxy composite plate specimen is mod- 
elled as a homogeneous transversely isotropic continuum plate medium. 
Acousto-ultrasonic non-contact input-output characterization by trac- 
ing SH waves in the continuum is studied theoretically with a trans- 
mitting and a receiving transducer located on the same face of the 
plate . 

The single reflection problem at a stress -free plane boundary in a 
semi- inf inite transversely isotropic medium whose isotropic plane is 
parallel to the plane boundary is analyzed first. It is found that an 
incident SH wave results in a reflected SH wave only; the amplitude 
ratio of the reflected SH wave to the incident SH wave is negative 
one; and the angle of reflection of the reflected SH wave is equal to 
the angle of incidence of the incident SH wave. The balance in energy 
flux normal to the plane boundary is also checked. 

The delay time is calculated as if the SH waves were propagating 
in an infinite half space. It is found that the directional depen- 
dence of the phase velocity of the SH waves travelling in the trans- 
versely isotropic medium has a significant effect on the delay time, 
as opposed to the directional independence of the phase velocity of 
the SH wave travelling in an isotropic medium. 
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The displacement associated with the SH wave in the plate detected 
by the non- contact receiving transducer is approximated by an asymp- 
totic solution for an infinite transversely isotropic medium subjected 
to a harmonic point load. The polar diagrams for the directivity 
function of the shear stress due to SH waves in the plate are shown at 
frequencies of 0.75, 1.50 and 2.25 MHz. To maximize signal detection, 
the transverse transmitting transducer should be placed such that the 
applied force generated by it is in the direction perpendicular to the 
line connecting the transmitting and the receiving transducers. This 
study enhances the quantitative understanding of acousto-ultrasonic 
nondestructive evaluation (NDE) parameters such as the stress wave 
factor (SWF) and wave propagation in fiber reinforced composites or 
any other materials which can be modelled as transversely isotropic 
media. 
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INTRODUCTION 


Fiber reinforced composite materials are attractive materials 
for aerospace applications because of their high specific mechanical 
properties. It has been shown that many composites, such as the uni- 
directional fiberglass epoxy composites or fiber reinforced ceramics, 
as shown in Fig. 1, may be modelled as a homogeneous transversely iso- 
tropic continuum [1] . In this work, acousto-ultrasonic (AU) non- 
contact input-output characterization of a homogeneous transversely 
isotropic elastic plate is investigated by tracing SH waves. 

First, single reflection problem of an incident SH wave at a 
stress-free plane boundary in a semi-infinite transversely isotropic 
medium whose isotropic plane is parallel to the plane boundary is con- 
sidered. At such boundaries, the conditions for the existence of wave 
mode conversion, the angle of reflection of the reflected wave, and 
the amplitude ratio of the reflected wave to the incident wave are 
derived. 

Second, the SH wave input- output relations are derived when mul- 
tiple reflections occur at the top and bottom faces of the plate. The 
delay time between input and output versus the distance separating the 
transmitting and receiving transducers is analyzed. The directivity 
function of the shear stress associated with the SH waves is computed. 
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And, the output displacement at the non-contact receiving transducer 
is approximated by an asymptotic solution. 

This investigation should enhance the quantitative understanding 
of AU NDE parameter such as the stress wave factor. It also provides 
the potential for assisting in the development of better NDE schemes 
utilizing the SWF. 


A 


SINGLE REFLECTION PROBLEM AT A STRESS -FREE PLANE BOUNDARY IN SEMI- 
INFINITE TRANSVERSELY ISOTROPIC MEDIUM WHOSE ISOTROPIC PLANE IS PARAL- 
LEL TO PLANE BOUNDARY FOR INCIDENT SH WAVE 

1 . Reflected SH Wave 

For a homogeneous linearly elastic transversely isotropic con- 
tinuum, the number of independent elastic constants is five [1] . 

Define a coordinate system (x, y, z) for a semi-infinite transversely 
isotropic medium whose isotropic plane is parallel to the plane bound- 
ary where the reflection occurs as follows: the plane boundary con- 

tains the x and y axes , and the z axis is the zonal axis of the 
medium, which is in the direction parallel to the fiber direction shown 
in Fig. 1, as shown in Fig. 2. The generalized Hooke's law is writ- 
ten, relative to the (x, y, z) coordinate system, as [1]: 

r xx “ c ll u - x + c 12 v >y + c 13 w > z 

Tyy - C 12 U,X + C U V,y + CjjW.Z 

r zz “ c 13 u - x + c 13 v -y + c 33 w > z 

r xz “ c 44< u > z + w,x) (1) 

T yz “ c 44< v - z + w.y) 

r xy “ c 66< u -y + v * x > 
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where r rs (r, s - x, y and z) are the normal (r - s) and shear (r * s) 
stresses with respect to the coordinate system (x, y, z) ; u, v and w 
are the displacement components of a point in the medium along the x, 
y and z axes, respectively; " , " denotes partial differentiation with 
respect to the variable which follows; and C^, ^ 12 » ^13* C33, and C44 
are the five independent elastic constants where Cgg - 1 / 2 ( 0 ^ - C^) . 

Let a plane progressive wave be represented as [ 2 ] 

(u, v, w) - A (P x , P y , P z ) exp{iw(S x x + S y y + S z z - t) } ( 2 ) 

where S x , S y and S z are the components of the slowness vector, which 
is in the same direction as the normal to the wavefront and whose mag- 
nitude is equal to the reciprocal of the magnitude of the phase velo- 
city [ 1 ], along the x,y and z axes, respectively; P x , P y and P z are 
the components of a unit vector of particle displacement along the x, 
y, and z axes, respectively; A is the amplitude of particle displace- 
ment; t denotes time and w denotes radian frequency. It follows from 
Eqs . ( 1 ) and ( 2 ) that the stresses can be represented as 
r xx - iwA[C^^S x P x + Ci2 s yJ?y + C 1 3S Z P Z ] exp{ iw(S x x + S y y + S z z - t) } 

Tyy - iu)A[C 12 S x P x + C ll S y P y + C 13 S z P z ]exp{iw(S x x + S y y + S z z - t) } 

r zz - iw A[C 1 3S x P x + C 1 3S y P y + C3 3 S Z P Z ] exp( iu>(S x x + S y y + S z z - t) } 

r xz “ ia) M c 44 S z p x + C44S x P z ]exp(iw(S x x + S y y + S z z - t) } ( 3 ) 

T yz “ io>A[C44S z Py + C44S y P z ] exp { iw ( S x x + Syy + S z z - t) ) 

r xy ** ^ w A[CggS y P x + CggS x P y ] exp ( icj(S x x + S y y + S z z - t) ) 

The stress boundary conditions on the stress -free plane boundary 
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require that [2] 


r xz 

r yz 

r zz 


(I) 

(I) 

(I) 


+ - 0 

+ r yz ( R > - 0 


( 4 ) 


where r rz ^^ (r-x, y and z) represents stresses on the plane boundary 
associated with the incident wave, and r rz ( R ) (r - x, y and z) 
represents stresses on the plane boundary associated with the 
reflected waves. In order to satisfy Eq. (4), it is required [2] that 
the frequency, gj, of the reflected waves be equal to that of the 
incident wave and that 


s x a) - s x w 

(5) 

S y <» - S y < R > . 

As a result of Eq. (5), the slowness vectors of the incident and 
reflected waves lie in a plane called the plane of incidence. This 
analysis can be simplified by assuming that the slowness vectors of 
the incident and reflected waves are in the plane x - 0, as shown in 
Fig. 2. Since the slowness vectors of the incident and reflected 
waves are in the plane x - 0, 

s x (I) “ S X < R) - 0- (6) 
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It has been shown that an SH wave with slowness vector in a 
plane containing the zonal axis of a transversely isotropic medium 
possesses a transverse displacement; that is, for the coordinates in 
Fig. 2, (P x , P y , P z ) - (1, 0, 0) [3]. It follows from Eqs. (3) and 
(6) that the stresses associated with SH waves are 

r xx * r yy ” r zz “ r yz “ ® 
r xz ** ® ’ r xy * 

It has also been shown [3] that an SV wave and a P wave with slowness 
vector in a plane containing the zonal axis of the transversely iso- 
tropic medium are quasi- transverse and quasi-longitudinal, respec- 
tively; the component of the unit vector of particle displacement for 
SV and P waves along the x axis, P x , vanishes, whereas the components 
along the y and z axes, Py and P z , do not. It follows from Eqs. (3) 
and (6) that stresses associated with P and SV waves are 


T yz * 0 


r xy “ T xz 
r xx ** 0 * 



r zz * 


0 . 


( 8 ) 


Assume that an SH wave is incident on the plane boundary, the 
x-y plane in Fig. 2. It follows from Eqs, (4) and (7) that 
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Since a reflected SH wave results in nonzero values of shear 
stresses, r xz and r x y, it is therefore concluded that only an SH wave 
will be reflected back into the medium. 


2 . Slowness Surface for SH Wave 

The equations of motion relative to the coordinate system 
(x,y,z) are [2] 

r xx, x + r xy,y + r xz,z “ pu,tt 

r xy, x + r yy,y + r yz,z “ P v » tc (11) 

r xz, x + r yz,y + r zz,z " pw,tt 
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It follows from Eqs. (1), (2) and (11) that the following equa- 


tions of motion are obtained: 

t^ll^x^ + + C 44 S Z 2 ’ + (^12 + ^66^x^y^y 

+ (C 13 + c 44 )s x s z p z - 0 

(^12 + Cgg)S x SyP x + t^66^x^ + ^ll^y^ + ^44^z^ " ^]^y (12) 

+ (C 13 + C 44 )S y S z P z - 0 

(0^3 + ^44) ^x^z^x "*" (^13 + ^44) SyS z Py + [C 44 (S X 2 + Sy 2 ) 

+ c 33 s z 2 '^l p z ” 0 

The condition of the plane wave solution is expressed by setting 
the determinant of the matrix of the coefficients of P x , P y and P z in 
Eq. (12) equal to zero [1]: 


[CnS x 2 ♦ c^Sy 2 ♦ c 44 s z 2 -p] 

< C 12 * C 66> S x S y 

< c 13 * C 44 )S x S z 

CC 12 ♦ C^Sy 

tC66 s x 2 + C 11 S y 2 + c 44 s z 2 

(C 13 * C 44 )S y S * 

(C 13 + C 44 )S X S Z 

(C13 ♦ C44 )S y S z 

tC 4 4« x 2 * S y 2 ) ♦ C 33 S Z 2 -p] 



(13) 


Expanding Eq. (13), we obtain three sheets of slowness surface. The 
slowness surface for an SH wave is [3] 

C 66(S X 2 + S y 2 ) + C 44 S Z 2 - p . (14) 
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In the present study, the slowness vectors of the incident and 
reflected SH waves are confined to the y-z plane, as shown in Fig. 2. 
The intersection of the slowness surface with the plane x-0 is given 
by [3] as 



The numerical values of the constants C 44 , Cgg and p for the 
fiberglass epoxy composite shown in Fig. 1 are given in [1] as 
C 44 - 4.422 X 10 9 N/m 2 , C 66 - 3.243 x 10 9 N/m 2 , and p - 1850 kg/m 3 . 

One quadrant of the slowness surface for an SH wave travelling in the 
fiberglass epoxy composite is obtained by substituting the numerical 
values of the elastic constants and the density into Eq. (15), and is 
shown in Fig. 3. 

3 . Angle of Reflection 

It follows from Eq. (5) that the component of the slowness vec- 
tor of an incident SH wave along the y axis, Sy^^, is equal to that 
of the reflected SH wave, Sy( R ^ . As result of Eqs . (5) and (15), the 
relationship between the z- components of the slowness vectors of an 
incident SH wave and the reflected SH wave is 
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s z m . . Sz w. 


(16) 


The minus sign is due to the fact that the slowness vector of an 
incident SH wave points out of the medium, whereas the slowness vector 
of the reflected SH wave points into the medium, as shown in Fig. 2. 
Consequently, the value of the z-component of the slowness vector of 
an incident SH wave, S z ^\ is negative, whereas that of the reflected 
wave, S Z ( R >, is positive. 

The angle of reflection is defined as the angle between the 
slowness vector of a reflected SH wave and the normal to the plane 
boundary where the reflection occurs. Similarly, the angle of inci- 
dence is defined as the angle between the slowness vector of an inci- 
dent SH wave and the normal to the plane boundary, as shown in 
Fig. 2. Therefore, the angle of reflection of a reflected SH wave is 

$ R - tan- 1 (S y ( R )/S z < R >) (17) 

and the angle of incidence of an incident SH wave is 

9 r - tan- 1 (S y ( I >/-S z ( I >). (18) 

It follows from Eqs. (5), (16), (17) and (18) that the angle of inci- 
dence is equal to the angle of reflection, as shown in Fig. 2. 
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4 . Amplitude Ratio of Reflected Wave to Incident Wave 


The amplitude ratio of the reflected SH wave to the incident SH 
wave can be computed from the boundary condition oh the stresses, 

Eq. (4). Assume an SH wave of unit amplitude travelling in the plane 
x * 0 is incident on the plane boundary at the origin, as shown in 
Fig. 2. The shear stress, r xz ^ , associated with the incident SH 
wave on the plane boundary can be obtained from Eqs . (3) and (6) as 

r X z (I) - i W C 44 S z ( I ) . (19) 

Similarly, the shear stress, r xz ( R ), associated with the 
reflected SH wave on the plane boundary is 

r xz (R) “ iwA( R >C 44 S z W. (20) 

It follows from Eqs. (10), (16), (19) and (20) that the amplitude of 
the reflected SH wave is 

A( R > - -1. (21) 

Therefore, the amplitude ratio of the reflected SH wave to the inci- 
dent SH wave is negative one. 
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5 . Balance in Energy Flux Normal to Plane Boundary 


The balance in energy flux normal to the plane boundary must be 
satisfied and is expressed as [3] 

F 2 a> + F z ( r ) - 0 (22) 

where F 2 '*) and F Z ^ R ) are the components of the energy fluxes of the 
incident and the reflected SH waves along the z axis, respectively. 

The component of the energy flux associated with an SH wave 
along the z axis is [4] 

F z - C 44 aVs z . (23) 

It follows from Eqs. (16), (21) and (23) that the sum of the z-compo- 
nent of the energy flux of the incident SH wave and that of the re- 
flected SH wave is equal to zero. Accordingly, the requirement for 
the balance in energy flux perpendicular to the plane boundary where a 
reflection occurs, Eq. (22), is satisfied. 

The z -components of the energy flux of an incident SH wave of 
unit amplitude and the reflected SH wave, Eq. (23), in the semi- 
infinite transversely isotropic medium are shown in Fig. 4, exclusive 
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of the radian frequency w, as functions of the angle of incidence with 
the numerical value of the elastic constant C 44 given in [1] for the 
unidirectional fiberglass epoxy composite shown in Fig.l as C 44 - 
4.422 x 10^ N/m^ . The balance in energy flux normal to the plane 
boundary, Eq. (22), is also shown in Fig. 4, which is independent of 
frequency. 
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ACOUSTO -ULTRASONIC NON -CONTACT INPUT-OUTPUT CHARACTERIZATION OF FIBERGLASS 
EPOXY COMPOSITE PLATE 

It has been shown [1] that the unidirectional fiberglass epoxy 
composite shown in Fig. 1 may be modelled as a homogeneous trans- 
versely isotropic continuum. For the axes shown in Fig. 1, the iso- 
tropic plane of equivalent continuum lies in the midplane of the plate 
[1]. A cartesian coordinate system (x, y, z) is chosen so that the x-y 
plane is the isotropy plane; thus, the upper and lower surfaces of the 
plate are at z- h /2 and z- -h/ 2 , respectively, where h is the plate 
thickness. The properties of the equivalent continuum model of the 
unidirectional fiberglass epoxy composite plate are [ 1 ] 

h - 0.1 m 

c n - 10.581 x 10 9 N/m 2 

C 12 “ ^-679 x 10 9 N/m 2 

C 13 - 40.741 x 10 9 N/m 2 (24) 

C 33 - 4.422 x 10 9 N/m 2 

C 44 - 3.243 x 10 9 N/m 2 

p - 1850 kg/m 2 

Non-contact transmitting and receiving transducers are located on 
the same face of a fiberglass epoxy composite plate specimen, as shown 
in Fig. 5. The unidirectional fiberglass epoxy composite plate spe- 
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cimen is considered as a plate of thickness h and of infinite planar 
(x-y) extent. The input electrical voltage to the transmitting trans- 
ducer is V^(t) and the output electrical voltage from the receiving 
transducer is V Q (t) where t represents time. The transmitting trans- 
ducer converts an input electrical voltage into a stress, whereas the 
receiving transducer converts a stress associated with stress waves 
travelling in the plate into an output voltage. In the following ana- 
lysis, only the SH waves are traced. The SH waves which are generated 
by the transmitting transducer located above point 0 experience mul- 
tiple reflections at each face of the plate, and then reach the 
receiving transducer located above point M, as shown in Fig. 6. Since 
the isotropic plane lies in the midplane and is parallel to both the 
top and the bottom faces where multiple reflections occur, the angle 
of reflection of the reflected SH wave is equal to the angle of inci- 
dence of an incident SH wave for each reflection at each face of the 
plate. Accordingly, the SH wave travelling from the input 0 to the 
output M may be considered as waves propagating in a semi -infinite 
transversely isotropic medium and travelling to point M' as if there 
were no bottom face, as shown in Fig. 6. 

1. Delay Time and Phase Velocity 

Let the input 0 and the output M lie in the y-z plane. Assume 
the number of reflections at the bottom face experienced by the SH 
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wave in travelling from the input 0 to the output M is n, as shown in 
Fig. 6. With respect to the z axis, the angle of incidence of the SH 
wave at each face of the plate is 8 , and the total distance travelled 
by the wave is R^. From the geometry in Fig. 6, 

8 - tan _1 (i/2nh) (25) 

where Z is the separation distance between the input 0 and the output 
M, and 


R n - Z/ sin*. (26) 

The delay time t n for the wave to reach the receiving transducer is 

t n “ *n/Cl(*) (27) 

where C^(0) is the directionally dependent phase velocity of the SH 
wave, and which for the unidirectional fiberglass epoxy composite 
plate is [1] 

Ci(0) - [(C 66 sin 2 0 + C^cos 2 *] 1 ^ (28) 

where C44, Cgg and p are given by Eq. (24). 
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The phase velocity as a function of the angle of incidence 9 
is shown in Fig. 7. The delay time is then computed when the number of 
reflections n at the bottom face of the plate is equal to 10, 100, 300 
or 500. The numerical results are shown in Fig. 8 where the delay time 
t n is plotted as the ordinate, and the dimensionless quantity Z/h is 
plotted as the abscissa for values of zero to 300. 

2 . Displacement Detected Bv the Receiving Transducer 


The displacement detected by the non- contact receiving trans- 
ducer above point M, radiated by the transverse transmitting trans- 
ducer, is assumed to be equivalent to the displacement at point M' , 
associated with the SH wave propagating in a semi- inf inite trans- 
versely isotropic medium. Except for the reflection coefficients at 
each face (to be discussed later), the displacement is computed as if 
there were no bottom boundary, as shown in Fig. 6. The shear stress 
at point M' is approximated by the far- field asymptotic solution for 
large of an infinite transversely isotropic medium subjected to a 
harmonic point load. 

Now consider an infinite transversely isotropic medium in which 
the z axis of a rectangular cartesian system 0(x,y,z) is the zonal 
axis of the medium and the x-y plane coincides with the isotropic 
plane, as shown in Fig. 9. The equations of motion including the 
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body force are [5] 


r xx,x + r xy,y + r xz,z + pX - pu,tt (29) 

r xy ,x + r yy ,y + r yz,z + P Y “ pv.tt (30) 

r xz,x + r yz,y + r zz,z + pZ — pw,tt (31) 

where r rs (r, s « x, y, z) are the normal (r - s) and shear (r ^ s) 

stresses with respect to^the chosen coordinate system 0(x, y, z); u, 
v and w are displacement components of a point in the medium along the 
x, y and z axes, respectively; X, Y and Z are the components of the 
body force along the x, y and z axes respectively; p is the density; t 
is time; and " , " denotes partial differentiation with respect to the 
variable which follows. 

If we differentiate Eq. (30) with respect to x and Eq. (29) with 
respect to y, we find upon substitution of the resulting equations and 
using the appropriate stress-strain relations, Eq. (1), that 

A , tt - — (A, xx + A ,yy) + C |-- A,zz + Y,x - X,y (32) 

where A - v,x - u,y. 
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For a harmonic point load at the origin, we take body forces of 


the form 


X - X 0 5(x)S(y)*(z)e' iut 
Y - Y q 5 (x)S (y)S (z)e' ^ ut 
Z - Z 0 S(x)S(y)S(z)e- iut 


(33) 


where 6 (r) (r - x, y, z) are the Dirac delta functions and x 0 , Y q and 
Z Q are the magnitude of the respective point body forces. A can be 
expressed as a threefold Fourier integral [5] 


A(x,y , z , t) - 


A(S x ,Sy,S z ,t) exp{ iw(S x x+Syy+S 2 z-t) )dS x dSydS z (34) 


where 


A(S x ,S y ,S z ,t) - 1/8 w 


A(x,y,z,t) exp(i<j(-S x x-Syy-S z z+t) }dxdydz. 


Similarly, Y,x and X,y can be expressed as threefold Fourier integrals 


Y,x - iwS, 


Y exp( iw(S x x+Syy+S z z- t) )dS x dSydS z 


(35) 


X,y - iwS, 


X exp{iw(S x x+Syy+S z z-t) )dS x dSydS z 
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where 


X - 1/8*- 


X 0 5(x)S(y)5(z)e 


■ iwt 


exp{iw(S x x+Syy+S z z-t) JdS x dSydS z 


- X 0 /8* 3 


and, similarly. 


Y - Y 0 /8* 3 ; Z - Z 0 /8* 3 . 


It follows from Eqs . (32), (33), (34) and (35) that 
i(S x Y 0 - S X Q ) 

A (36) 

8* 3 wG(S x , Sy , S z ) 

where G(S X , Sy, S z ) ■* Cgg/p(S x ^ + Sy^) + C^^/pS z -1. 

Therefore, G(S X , Sy, S z ) - 0 represents the slowness surface for 
an SH wave by comparison with Eq. (14). As a result of Eqs. (34) and 
(36) , A becomes 


A(x,y,z,t) - 


iSxYo-iSyXo 

8x 3 wG(S x ,Sy,S z ) 


exp(iw(S x x+Syy+S z z-t) JdS x dSydS z (37) 
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The asymptotic solution at large distance from the point load is 
obtained by applying the theory of residues, the method of stationary 
phase, and the radiation condition [5] as 


iw(S *Y_ - S *X„) * * * 

A(x,y,z,t)-. X- exp(iw(S*x+S *y+S z *z-t)) (38) 

2ttR|vg| ( | Kq | 

^ 

where (S x , Sy , S z ) are the points on the slowness surface for an SH 
wave, G (S x , Sy, S z ) - 0, where the normal is parallel to the OP 
direction; R is the distance from the origin 0 where the point load is 
applied to the point of interest P in the medium; VG is the gradient 
of G (S x , Sy, S z ) and | VG | is the magnitude of the gradient of G(S X , 
Sy, S z ) and is expressed as 


|VG| - [G,s x + G,Sy + G,S z ] 1 /2 


(39) 


and is evaluated at points (S x *, Sy*, S 2 *) of the slowness surface for 
an SH wave, G(S X , Sy, S z ) - 0, where the normal is parallel to the OP 
direction; Kq is the gaussian curvature of the slowness surface for an 
SH wave and is expressed as 
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1 


k g 


(G?S X + G?S + G? S z ) 


S{ G , (G , S X S X G , SySy-G , S X Sy) 


+ 2G,S X G,Sy(G,S x S z G,S y S z -G,S x Sy G,S Z S Z )} 


(40) 


where 2 is the sum with respect to cyclic permutation of S x , Sy and 
S z> and ( | K g | ) 1/2 is the square root of the magnitude of the gauss ian 
curvature of the slowness surface, G( S x , Sy, S z ) - 0, and is evaluated 
at points (S x *, Sy*, S z *) of the slowness surface where the normal is 
parallel to the OP direction, as shown in Fig. 9. 


The displacement component along the x and y axes, u and v, due 
to an SH wave can be obtained by direct integration of the definition 
of A in Eq. (32), in conjunction with pages 14 and 15 of [6] and is 
asymptotically expressed as 


u(x,y,z, t) 


v(x,y,z, t) 


q *2y _ q *q ^y 

a y A o ^x D y 1 o 
2*R(S X * 2 + S y * 2 )|VG|(|K G |) 1 /2 

★ * * 

exp(iu>(S x x + S y y + S z z - t) 

-S x *Sy*Xo + S x * 2 Y q 
2ttR(S x * 2 + S y * 2 )|VG|(|K G |) 1 /2 
* * * 

exp{iw(S x x + Syy + S z z - t) 


(41) 
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It is known from the definition of A in Eq. (32) and Eq. (37) 
that SH waves have no contribution to the displacement component along 
the z axis, w. 


However, it should be noted that there is only one point on the 
slowness surface, G(S X , Sy, S z ) - 0, where the normal is parallel to a 
given direction OP. This is because the slowness surface for an SH 
wave is a spheroid, Eq . (14). Consequently, if (x, y, z) are the 
coordinates of a given point in the medium, there will be only one SH 
wavefront passing through it corresponding to a point (S x , Sy , S z ~) 
on the slowness surface where the normal is parallel to the given 
direction. Note that in this analysis, it has not been necessary to 
use Eq. (31) . 


3 . Directivity Function 


The shear stress, r xz , associated with SH waves reaching the 
point M' in Fig. 6 is given from Eqs . (1) and (41) as 


r xz (x,y,z,t)~ 


iC 44W S z *(S y * 2 X 0 - S x *S y *y o ) 
2;rR(S x * 2 + S y * 2 )|VG|(|K G |) 1 /2 


exp{iw(S x *x + S y *y + S z *z - t) } (42) 
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The directivity function D SH of the shear stress, r X2 , 
associated with the SH wave is determined from the amplitude of the 
shear stress in Eq. (42) by setting X Q - Y 0 - R - 1; 


C44S Z *|S y * 2 • S x *Sy*[q, 

2*(S X * 2 + Sy* 2 ) | | ( | Kq | ) L / 2 


(43) 


Since the point 0 and the point M' in Fig. 6 are located in y-z 
plane, the evaluation of the directivity function Dgfl is simplified by 
setting S x - 0 in Eq. (43). 


Due to axial symmetry with respect to the zonal axis, the z 
axis, of the transversely isotropic medium, the value of the directi- 
vity function thus obtained holds for all values of S x . 

Therefore, it follows from Eq. (39), (40) and (43) that 


d SH “ 


C44wS 2 


4 * c 66/^ c 44/*>) 1>/2 L f> 


c 66 * 0 c 44 * « 

( S y *) 2 + ( S z *) 2 


1/2 


(44) 


The polar diagrams for the directivity function of the shear 
stress, r xz , associated with the SH waves travelling in the unidirec- 
tional fiberglass epoxy composite are obtained from Eqs. (24), (25), 
and (44) at frequencies of 0.75, 1.50 and 2.25 MHz, and are shown in. 
Fig. 10, where the angle of incidence 0, given by Eq. (25), is used to 
determine the direction in which the value of a point (0, Sy*, S z *) on 
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the slowness surface where the normal is parallel to the given direc- 
tion is thus obtained. 

4 . Assumptions on the Transducers 

The non-contact transmitting transducer in Fig. 5 is assumed to 
transform an electrical voltage into a uniform shear stress; however, the 
non- contact receiving transducer in Fig. 5 transforms a displacement into 
an electrical voltage. The approach which follows is similar to that 
given in [7]. 

If an input voltage of amplitude V and frequency w is applied 
according to 

V L (t) - Ve’ iwt (45) 

see Fig. 5, the stress a that is introduced into the specimen by the 
non- contact transmitting transducer is 

o(t) - F 1 (w)Ve* i(wt + (46) 

where F^(w) is the transduction ratio for the transmitting transducer 
in transforming a voltage to a stress and ^ is a phase angle. In 
Eqs . (45) and (46), the harmonic character of the signals is expressed 


in the complex notation where i- J- 1 and only the real parts of these 
and subsequent equations should be considered. Thus, the amplitude T 
of the applied force is defined as 

T - F 1 (w)V. (47) 

Similarly, if a stress wave producing a displacement component U 
and frequency u> that is detected by the non- contact receiving trans- 
ducer is defined as 

u(t) - Ue iwt (48) 

the output voltage from the non-contact receiving transducer, see 
Fig. 5, is 

V Q (t) - F 2 (w)U' e* i ^ wt + (49) 

where F 2 (w) is the transduction ratio for the non- contact receiving 
transducer in transforming a displacemerft to a voltage, and <f >2 is a 
phase angle. Thus, the amplitude V' of the output electrical voltage 
is 


V' - F 2 (w)U. 


(50) 
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The characteristics of F^(w) and F 2 (w) are unknown except that the 
dimension of the product F^(w)F 2 (t*>) is [kg/m^ • sec^] . 

5 . Steady-State Output Voltage Amplitude due to Multiple Wave Reflec- 
tions in Plate 

Since only the SH waves travelling in the y-z plane are consid- 
ered in Fig. 6 , it follows from Eqs. (26), (39), (40) and (41) that 
the amplitude of the hypothetical displacement at point M' is U M ' , and 
is defined as 


(C 44 P ) 1 / 2 S z * T 

U M , (51) 

2jt Rn (C 11 -C 12 ) 

where T represents the magnitude of the shear stress generated by the 
non- contact transmitting transducer along the direction perpendicular 
to the line connecting the point 0 and the point M in Fig. 6 and is 
equivalent to X Q in Eq. (41). 

However, with the bottom boundary present, the wave is reflected 
a total of (2n-l) times, as shown in Fig. 6 . Thus, the amplitude of 
the displacement at point M is U^, and is obtained by modifying 
Eq. (51) as 


2 n-l 

U M “ Qshsh U M' 


(52) 
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where Qshsh *- s re f lection coefficient of SH waves to SH waves, and 

its value is -1, as given by Eq. (21), 

The amplitude of the output voltage from the non-contact receiv- 
ing transducer is V' and can be obtained by combining Eqs. (50), 
(51)and (52) as 


F 2 (*> QSHSH (C44P) 1/2 S z * T 

V' (53) 

2 * Rn ( c ll- c 12> 


Substitution of Eq. (47) into Eq. (53) gives 


2n-l 

F L (w) F 2 (w)Q shsh (C 44 p)V2 s z * V 

V' . (54) 

2n Rn (C 11 -C 12 ) 


Introducing the SH wave attenuation constant a of the fiberglass 
epoxy composite and a possible electrical signal amplification factor 
K, Eq. (54) can be written as 


K FiM F 2 (u) Q SHSH (C 44 p) 1 /2 s z * v e- aRn 

V' . (55) 

2k Rn (C^-C^) 


Eq. (55) gives the output voltage amplitude from the non-contact 
receiving transducer due to an input voltage amplitude V at the non- 
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contact transmitting transducer when the SH wave path has included n 
reflections from the bottom face of the fiberglass epoxy composite 
plate specimen, as shown in Fig. 6. 
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DISCUSSION AND CONCLUSION 


In the acousto-ultrasonic input-output characterization of the 
unidirectional fiberglass epoxy composite plate, the angle of reflec- 
tion of the reflected SH wave is equal to the angle of incidence of 
the incident SH wave for each reflection at either the top or the bot- 
tom face of the plate. This is due to the fact that the isotropic 
plane is parallel to either face of the plate. However, if the paral- 
lelism between the isotropic plane and the plane boundaries where 
reflection occurs does not exist, the angle of reflection is not equal 
to the angle of incidence. In such a case, the use of a semi- inf inite 
transversely isotropic medium, neglecting the existence of the bottom 
face of the plate (except for the cumulative reflection coefficients) , 
to compute the delay time t n and the directivity function becomes 
inappropriate . 

It is observed from Fig. 8 that an increase in the number of 
reflections n from the bottom face of the plate results in a minor 
increase in the delay time t n at each value of dimensionless separa- 
tion i/h. For a given plate thickness h and a given separation dis- 
tance £ between the transmitting transducer and the receiving trans- 
ducer, an increase in the number of reflections n results in a 
decrease in the angle of incidence 0, as given by Eq. (25), thereby 
increasing the travelling distance 1^, from Eq. (26). The phase 
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velocity of an SH wave in the unidirectional fiberglass epoxy 
composite increases with decreasing angle of incidence 6 as shown in 
Fig. 7. The increase in the travelling distance is partly offset 
by the increase in the phase velocity with an increase in the 
number of reflections n. This accounts for the minor increase in the 
delay time given by Eq. (27). On the other hand, in an isotropic 
medium, the phase velocity of a stress wave is directionally indepen- 
dent. Thus, the increase in the delay time t n caused by an increase 
in the number of reflections n is solely attributed to the increase in 
the travelling distance in an isotropic plate. Therefore, it is 
concluded that the directional dependence of the phase velocity of a 
stress wave travelling in a transversely isotropic medium has a sig- 
nificant effect on the delay time when conducting acousto-ultrasonic 
testing. 

Now, consider the case of a given number of reflections n. An 
increase in the separation distance i results in an increase in the 
travelling distance 1^, as shown in Fig. 6, and an increase in the 
angle of incidence $ given by Eq. (25). The increasing angle of 
incidence 6 results in the decreasing phase velocity C^ f as shown in 
Fig. 7. Consequently, the delay time t n increases sharply as a result 
of the increasing separation distance i for a given number of reflec- 
tions n, as shown in Fig. 8. However, as the separation distance Z 
approaches infinity , the phase velocity approaches a limit. Thus, 
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the increase in the delay time t n is due solely to the increase in the 
travelling distance R^, which is similar to the isotropic medium case. 

An examination of Eqs . (41) shows that the y-component displace- 
ment v associated with an SH wave travelling in the y-z plane, S x - 0, 
is zero and the x- component displacement associated with an SH wave 
travelling in the y-z plane, S x - 0, is attributed solely to the 
x-component of the applied point load. In other words, the non- 
contact receiving transducer is unable to detect the displacement com- 
ponent along the line connecting the transmitting and the receiving 
transducers by tracing SH waves in the plate. Further, the only 
detectable displacement component which is perpendicular to the normal 
to both faces of the plate and to the line connecting the transmitting 
and the receiving transducers (that is, the u displacement) becomes 
non-detectable if the shear stress introduced by the non-contact 
transmitting transducer is in the direction parallel to the line con- 
necting the transmitting and the receiving transducers, regardless of 
how large the stress is. Accordingly, to maximize the displacement 
detected by the non- contact receiving transducer, the transmitting 
transducer should be placed such that the shear stress generated by it 
is in the direction perpendicular to the line connecting the transmit- 
ting and the receiving transducers as shown in Fig. 11. 
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This theoretical investigation provides a step forward in the 
quantitative understanding of acousto-ultrasonic non-destructive eval- 
uation (NDE) parameters such as the stress wave factor (SWF) in trans- 
versely isotropic media. It also provides the potential for assisting 
in the development of more efficient and more revealing NDE schemes 
utilizing wave propagation. 
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X 


Coordinate system (x,y,z) in analysis of single reflection 
problem at stress-free plane boundary of semi- inf inite 
transversely isotropic medium; z = 0 is plane boundary 
where single reflections occurs, z < 0 is free space, and 
x = 0 is plane of incidence. 
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Fig. 3 Slowness surface for SH wave in unidirectional fiberglass epoxy c 
positive y-z quadrant. 



Fig. 4 Energy fluxes normal to plane boundary, exclusive of radian frequency, of incident 
and reflected SH waves in unidirectional fiberglass epoxy composite. 




Schematic of non-contact acousto-ultrasonic test configuration 
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- . . . wave which arrives at point 
tions from bottom boundary. 
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Delay time versus dimensionless separation £/h with number of reflections n 
from bottom face of unidirectional fiberglass epoxy composite plate specimen 
as parameter. 




Schematic illustrating harmonic point load exciting 
an infinite transversely isotropic medium, where xy 
is isotropic plane in cartesian coordinate system 
defined by (x,y,z). 
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Fig. 10 Polar diagram of directivity function of shear stress T xz associated 

with SH waves at frequencies of 0.75, 1.50 and 2.25 MHz for positive y-z quadrant 



.Direstion in which shear stress generated 
by transmitting transducer wilt result in 
maximum displacement detected by the 
non- contact receiving transducer by 
tracing SH waves 



Receiving 

transducer 


Transm itting 
t ransducer 


Fig. 11 Schematic illustrating maximization of displacement 

detected by non-contact receiving transducer by placing 
transmitting transducer such that shear stress generated 
by transmitting transducer is in direction perpendicular 
to line connecting transmitting and receiving transducers. 
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